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We investigate the role of equilibrium and nonequilibrium noise in the magnetization dynamics 
on mono-domain ferromagnets. Starting from a microscopic model we present a detailed derivation 
of the spin shot noise correlator. We investigate the ramifications of the nonequilibrium noise on 
the spin torque dynamics, both in the steady state precessional regime and the spin switching 
regime. In the latter case we apply a generalized Fokker-Planck approach to spin switching, which 
models the switching by an Arrhenius law with an effective elevated temperature. We calculate the 
renormalization of the effective temperature due to spin shot noise and show that the nonequilibrium 
noise leads to the creation of cold and hot spot with respect to the noise intensity. 
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I. INTRODUCTION 

The manipulation of magnetization of a ferromagnet 
by means of spin-polarized currents is a key issue of 
the state-of-the-art spintronics concepts (for a review see 
Ref[l.]). In this respect the most important phenomenon 
is the so-called spin-transfer-torque (STT) effect, which 
was predicted by Slonczewski and Berger 2,3 : A spin- 
polarized current may transfer angular momentum to a 
free ferromagnetic layer resulting in a macroscopic torque 
on the latter's magnetization. In very small ferromag- 
nets, in which the magnetization can be assumed spa- 
tially uniform, the STT results in the rotation of the 
magnetization as a whole rather than in the excitation of 
spin waves. The STT can in particular lead to two dy- 
namical regimes: the reversal of the free layer's magneti- 
zation or a steady state precession of the magnetization. 
Due to the giant magnetoresistance effect, the dynamics 
of magnetization is reflected in the change of resistance 
of the circuit containing magnetic junctions. Current in- 
duced resistance variations, identified with STT, were re- 
ported in Refjj]. In subsequent years both spin torque 
induced magnetization reversal were observed^ as well 
as steady state precession 7 - - — . Both dynamical regimes 
are interesting for applications: as a tempting alterna- 
tive to Oersted fields in switching the magnetization in 
ferromagnetic memory elements on the one hand, and as 
clock devices used to synchronize the CPU with other 
logic units, on the other. Hence a reliable description of 
these phenomena is very important. 

On a semi-classical level, magnetization dynamics of a 
mono-domain ferromagnet can be well described by the 
Landau-Lifshitz-Gilbert Equation (LLG) 

dm dm 70 , T . 

— — = -70m x H cff + a m x — - + -7777 m * (m x I s ) ■ 
at at M S V 

(1) 
Here m is a unit vector in the free layer's magnetization 
direction, H e g the effective magnetic field, V the vol- 



ume of the switching element, M s the absolute value of 
the free layer's magnetization, 70 the gyromagnetic ra- 
tio, oiQ the Gilbert damping parameter, and I s the spin 
polarized current. However, as the extension of associ- 
ated devices are very small, effects of noise may play a 
significant role and should be included into the dynam- 
ical description. A first inclusion of noise into the LLG 
was given in the seminal paper of Browni 1 -, who consid- 
ered the effect of thermal fluctuations on the dynamics 
of a mono-domain particle by a random component h^ 
of the effective magnetic field entering the LLG equa- 
tion ([1]) . As a consequence of the fluctuation-dissipation 
theorem one finds for the equilibrium correlator of the 
random field 11 



(Jii(t)fcj(O) « a Q k B T5 l3 5{t - t') 



(2) 



where hi(t) denotes the i-th Cartesian component of the 
random field at time i, ks is the Boltzmann constant 
and T the temperature. Since then, temperature ef- 
fects on the LLG equation have been considered, both 
witbJ£~— and without the spin torque ter m 17 ' 18 . The 
emphasis of these approaches has been on the influence 
of noise on switching rates, often by performing explicit 
numerical calculation a 12 i 13 ' 17 i 18 . Since the spin-torque 



experiments 



8.9.19-21 



are performed under clear nonequi- 



librium conditions, it is natural to address, besides tem- 
perature, other sources of noise. One possible source of 
nonequilibrium noise is the spin shot noise. By analogy 
with the charge shot noise the quantization of the angular 
momentum transfer leads to spin shot noise, which mani- 
fests itself in a random Langevin force entering the equa- 
tions of motion for the free magnetic layer. It was shown 
by Foros et aW* in the context of normal metal / ferro- 
magnet / normal metal (NFN) structures that the spin 
shot noise is the dominant contribution to magnetization 
noise at low temperatures. In the realistic experiments 
on spin torque and spin switching the nonequilibrium 
noise starts to dominate at temperatures below several 
Kelvins. In Ref 221 it was shown that inclusion of the 



nonequilibrium noise into the dynamical description can 
explain the experimentally obseved nonmonotonic depen- 
dence of the microwave power spectrum on the voltage, 
as well as the saturation of the spectral linewidth at low 
temperatures. 

In this paper we concentrate on the effect of the 
nonequilibrium noise on spin-switching. Without noise, 
the spin switching takes place when the spin current 
(and with it the resulting spin torque) exceeds a critical 
value. The critical current is determined by the strength 
of the magnetic anisotropy, external magnetic field and 
Gilbert damping, and it can be obtained from the solu- 
tion of the deterministic LLG equation ([1]). Magnetiza- 
tion noise opens a possibility of activated switching at 
currents much less than critical. Moreover, in the pres- 
ence of noise the switching becomes a random process 
that requires a probabilistic description. The latter is 
based upon the solution of the Fokker-Planck (FP) equa- 
tion for the probability distribution of magnetization as 
derived in Ref[22|]. A crucial step towards this descrip- 
tion has been made in Ref [lj] , where the authors reduced 
the FP equation for magnetization to the effectively one- 
dimensional FP equation for the probability distribution 
of energies and applied it to the description of activated 
switching processes. 

A number of experiments on current induced switch- 
ing have been carried out previously 2 ^ - — . From them 
it was found that the noise reduces the typical switch- 
ing time, as one might expect. Myers et al<^ observed 
a broad distributions of switching currents strongly de- 
pending on temperature, indicating a thermally acti- 
vated switching process altered by the STT. To fit the 
measured data they used a Neel-Brown mode l 11 ' 28 with 
a field- and current-dependent potential barrier height 
U(H,I). In this model the probability for the magne- 
tization to switch decays exponentially with time over 
a characteristic relaxation time r that obeys the Ar- 
rhenius law r ~ e u/k B T _ ^ n i m pii c it assumption in 
the Neel-Brown theory is that magnetization dynamics 
is governed by a torque from an effective magnetic field 
H e g, which is derivable from the free energy E(M.) of 

the system via H e g = Vj\f£(M). The spin torque 

however is non-conservative and the concept of a cor- 
responding potential barrier is ill-defined, which makes 
the situation significantly more complicated. For ther- 
mally activated switching in presence of STT Urazhdin et 
a L 24 ^ 5 found that the activation energy strongly depends 
on the magnitude as well as the direction of the cur- 
rent. To capture the observed features they introduced 
an effective temperature distinct from the real tempera- 
ture in the Neel-Brown formula. Its current directional 
dependence indicated that the heating is not the ordi- 
nary Joule heating. Based on a stationary solution of 
the Fokker-Planck equation Apalkov and Vissche r 14 ' 15 , 
and Li and Zhang 12 , linked this effective temperature to 
the spin torque. In their model the alteration of switch- 
ing rates is due to the change of the elevated effective 
temperature in the Arrhenius factor, which in general 



yields a non-Boltzmann probability distribution. We ex- 
tend the approach of Ref [14j considering the influence of 
the noise on the switching in the whole range of currents, 
from the noise induced activated switching at small cur- 
rents up to the almost deterministic switching by large 
critical currents. We also take into account specific angu- 
lar dependence of the nonequilibrium noise and analyze 
the applicability of the effective temperature description 
to the nonequilibtium noise in detail. 

With the present paper we hope to give a contribution 
towards a better understanding of noise in magnetic sys- 
tems. In sectionlLTlwe start our discussion with a detailed 
derivation of the spin shot noise correlator by means of 
the Keldysh technique. It is our goal to depict the under- 
lying mechanisms that lead to the occurrence of spin shot 
noise in magnetic nanodevices and to provide a general 
mathematical framework for their description. We then 
turn our attention to the ramifications of the noise on 
the magnetization dynamics. We address the question 
of switching rates estimation in sections IIIII and IIVI by 
applying a generalized Fokker-Planck approach. Within 
this approach the alteration of switching rates due to 
spin torque is determined by an effective temperature 
T e g. The latter differs from the real temperature T, as it 
incorporates the effects of the damping, the spin torque, 
and noise. We calculate the renormalization of the ef- 
fective temperature due to the nonequilibrium noise. A 
conclusion of our findings is given in section [Vj 



II. SPIN SHOT NOISE CORRELATOR 

In this section, starting from a microscopic model of a 
magnetic tunnel junction (MTJ) we will derive a stochas- 
tic version of the LLG equation. Fluctuations will nat- 
urally come about due to the nonequilibrium situation, 
and will comprise the random part of the stochastic LLG. 
In particular, performing a perturbative expansion of the 
Keldysh action in terms of the spin flip processes and the 
tunneling amplitude we will be able to derive the spin 
shot noise correlator. 

The model MTJ consists of two itinerant ferromag- 
nets separated by a tunneling barrier. Let us introduce 
the corresponding model Hamiltonian allowing for an ex- 
ternal magnetic field H, tunneling of itinerant electrons 
through the barrier and exchange coupling between the 
itinerant electrons and the free layer's magnetization. It 
reads 



Ho = y^ efco-cj^Cfco. + 

fe.CT 



£%i4,*a-7S-H-2JS-s 

la 



J2 W ^ 



ka 



+ h.c. 



Ida 



(3) 



The notation is as follows: The creation (annihilation) 
operators c* ka (cfc CT ) and d\ a (d/ CT ) describe the itinerant 
electrons of the fixed and the free magnetic layer respec- 
tively, cr = + corresponds to the respective majority and 



a = — to the minority spin band, and the indices k and I 
label momentum. The operator S describes the total spin 
of the free layer. It is connected to the free layer's mag- 
netization via S = MV/7. s = i 'Yliiaa 1 d\ a o aa idi a i is the 
quantum operator associated with the spin of itinerant 
electrons, where a denotes the vector of Pauli matrices. J 
is the exchange coupling constant and Wki are tunneling 
matrix elements. 

For the subsequent discussion we assume that the time 
between two tunneling processes is much larger than 
the relaxation time in the free ferromagnet, which is 
equivalent to assuming a complete spin relaxation in the 
free magnetic layer. This allows us to introduce an in- 
stantaneous reference frame with the spin quantization 
axis directed along the free layer's magnetization direc- 
tion. To render the free layer's magnetization a dy- 
namical variable, we make use of the Holstein-Primakoff 
paramctrization 29 
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FIG. 1: (Color Online) Diagrams for spin flip processes: a) 
First order, b) Second order. Solid (dashed) lines denote 
electronic propagators in the free (fixed) layer. Bold dashed 
lines are propagators of HP bosons. Tunneling vertices are 
denoted by circles with crosses. 



the free and fixed layer we obtain in the energy domain 



S z 



S-tfb; S- =bW2S-tfb; S+ = ^2S- 



6t66, 

(4) 
where tf , b are usual bosonic operators and S± = S x ± 
iS y . At low temperatures we can assume that the expec- 
tation value of b'b is much smaller than 2S allowing to 
treat the square root to zeroth order in tfb. Taking all of 
the above mentioned into account, Hamiltonian Q can 
be written in the instantaneous reference frame as 



fe.cr 



^ka^ka 



(ei - JSa)d\ a di a - ^SH, 



la 



- 1 tfbH z + jtfbY, (jd \A« + 

la 

-b^2S[jY,d\ i d^ + 1 -H^ 



E wzrcidi*' 



r'J 



A/,, 



h.c. 



(5) 



where we used the notation H± = H x ± iH y . W£f are 
spin dependent tunneling matrix elements given by 



IF; 



ki 



{<j\a')W k i, 



(a\<r) = cos |e- = ff * 



[a\a ) — a sin^e 2 



(G) 

(7) 



Hamiltonian ([5]) can be now readily translated into a 
Keldysh action using the general scheme of the Keldysh 
technique 30 . To this end we switch to symmetric ("cl") 
and antisymmetric (" g" ) linear combinations of the field 
operators. In accordance with parametrization (U) the 
former are connected to the m± components of the 
free layer's magnetization in the instantaneous reference 
frame via 



qR/A 

la 



1 



e — eia- ± zO 



G 



R/A 

key 



1 



e - e ka ± iO 



. (9) 



where ti a = e; — a JS are the energies of the itinerant 
electrons with momentum I and spin a in the free ferro- 
magnet, and eka the corresponding energies for the fixed 
layer. The Keldysh components are 

Gf a = (l-2n d F (e))S(e-e la ), (10) 

Gl = (l-2n c F (e))5(e-e klT ), (11) 

where chemical potentials fJ-d/c f° r the free and fixed layer 
are included in the fcrmionic distribution functions 



c/d 



For future use we also define the matrices in Keldysh 
space 



7* = 



1 
1 



r 



1 

1 



(12) 



The Keldysh action can be now solved perturbatively 
in terms of the tunneling amplitude and the spin flip 
processes. In second order in both quantities this leads 
to the diagrams of Fig. [TJ The corresponding equations 
of motion are obtained when varying the action A with 
respect to the quantum component 



8A_ 5A 

Sb a Sb n 







(13) 



Finally we note that, in the instantaneous reference 
frame, we have for m± = m x ± im y 



(m±) - , 



(dtm±) £ 



(14) 



bd(t) = 



>M 8 V 



m+{t) , b c i(t) 



'M 8 V 



-(*). 



(8) 
For the retarded and advanced components of the 
fermionic Green functions for the itinerant electrons of 



where, on the other hand 



= 1, 



(d t m z ) = . 



(15) 



We may now translate diagrams Q] (a) and (b) into the 
analytical expressions. However, let us start with the 



contribution of zeroth order (in spin flips and in tunnel- 
ing). It reads 



dtb q (t) (id t b c i(t) + -f^S/2^ 



The resulting equations of motion are 

id t b c i+ 1 ^S]2H+=Q 



(16) 



(17) 



and a corresponding complex conjugate equation for b c i. 
Equation (fTT)) describes the precession of the magneti- 
zation around the magnetic field H and forms the first 
term of the LLG equation ([1]) . 

Let us come to the diagram of Fig. QJa). To extract its 
contribution to the action we have to calculate 

-JVSJ2 W^-'W^'baTriGt-fGf^G^} , (18) 

klcra' 

where for brevity the symbolic notation b a with fy- = b q 
and bi = b q was introduced. The resulting action reads 



A x 



2VS 



dt{b q (t) 



sin o e 



^ -b q (t) sin 9e irt> } 



(19) 

Variation of (|19p with respect to b q and b q gives the fol- 
lowing contribution to the equations of motion 



SAi 
Sb q (t) 



,i(t> 



25 



SAi 

5b q (t) 



25 



,-i<t> 



(20) 



Again, using the HP parametrization (J4J and the relation 
between S and m, equation (|2U1) can be readily translated 
into the corresponding equation of motion for the mag- 
netization. The result is the spin torque term of equation 



ftm 



7 



M S V 



m x (I s x m) 



(21) 



As far as the remaining diagram (Fig. |TJb)) is con- 
cerned we have to distinguish two contributions: One 
with two quantum components and one with a quantum 
and a classical component respectively. 36 In the first case 
we obtain 



J 2 Sb q b q J2 |W^-'| 2, &{Gf_ <7 ( e )7 9 Gt( e - W )x 

klcra' 

^GtMGlAe)} ■ (22) 
In the second case we have 

J 2 Sb cl b q £ \W^-"\ a Tr{Gf_ a (e)^Gt{ e - U )x 

7 d G?_ CT ( £ )G^(e)} , 

(23) 



klo 



and the corresponding contribution with q ■H- d. The 
resulting action is 



A 2 = 

where 
a(9) 



dt 



a(9) (b q d t b cl - bcidtbg) + ^V(9)b q b q 



% (dl.fffi) 



eMV V dV 



V(9) = ^ ao k B T+h sf (9) coth :T } 



(S=\ 



(24) 

(25) 
(26) 



The spin flip current I s f can be calculated from the elec- 
tric conductances GpiAP) in the parallel (antiparallel) 
configuration as follows 



dl sf (9) _ h 
dV 4e 



Gp sm~ I - 
1 2 



Gap'"'^ I - 



(27) 



Action (|2"4"|) consists of two parts. The first term is a 
damping term. In the LLG equation it will result in a 
renormalization of the Gilbert damping parameter. The 
renormalization is due to the coupling to the reservoirs. 
The enhancement of the damping, Eq. (|23|) . is closely 
related to the spin-pumping enhanced damping as dis- 
cussed in Ref 31.32| in the framework of the Landauer- 
Biittiker formalism. As far as the second term of Eq . (f2"4")l 
is concerned we introduce a Hubbard-Stratonovich auxil- 
iary field which decouples the action. Let us denote this 
(complex) field by if = lf x + ilfy We can write 






iAi 



dlfdlfe-^'+'+ e 

= j 'd/«d7«e-^+'+ exp{_i-^(J«6 B + J«6 g )} , 

(28) 

where we abbreviated the second term of ([24| by ^22- 
As one can see the result is a noise-averaged term which 
is linear in the quantum component. The linear action 
constitutes a resolution of functional (5-functions of the 
Langevin equations on b c i(t) and its complex conjugate. 
The stochastic properties are encoded in the auxiliary 
field /:?, precisely in the correlator ([26)) . For b c i the 
Langevin equations read id t b c i — -4=1^. This corre- 
sponds to idtm+ = jftfl+ leading to the random term 
of the stochastic LLG equation. Adopting the notation 
1^ = $!„, in conclusion we have found 



dtm 



7 



M S V 



m x (SI S x m) , 



(29) 



where the stochastic field is characterized by 

(8I a ,i(t)tt.j(fi)) = ^{0)SijS(t - t') (30) 

with the correlator T>(9) given by Eq. ([2"B)) . 

To complete our discussion we add some comments 
concerning the correlator ([2"B)) . To start with, we note 



that T> contains two parts, an equilibrium part (which 
is phenomcnological, and in compliance with the FDT 
proportional to ao taking into account intrinsic damping 
processes) and a nonequilibrium part. The nonequilib- 
rium part exhibits a dependence on the mutual orienta- 
tion of the fixed and free layer's magnetizations. This 
angle dependence enters the correlator through the spin 
flip current I s f . The physical meaning behind this quan- 
tity is the following: I s f counts the total number of spin 
flip events, irrespective of their direction. Hence, even 
if there is no contribution to the spin current J s , the 
spin flip current J s / may acquire a nonzero value. The 
discreteness of angular momentum transfer in each spin 
flip event leads to the occurrence of the nonequilibrium 
noise. In this sense the nonequilibrium part of (|26[) can 
be identified with the spin shot noise. 

In conclusion we have derived the following stochastic 
LLG equation 

dm xx dm 7 r , .du 

— = -7omxH cff +a mx— + — — mx m x (I s + If ) 
at at M S V 

(31) 



where the random field correlator is given by Eq. ([30 
along with ([H and (JS7) • 



III. FOKKER-PLANCK APPROACH TO SPIN 
TORQUE SWITCHING 

Spin torque switching is observable in two different 
regimes. On the one hand, the spin torque can switch 
the magnetization of a free ferromagnct when the current 
exceeds a critical value I c . On the other hand, switch- 
ing is also observed for currents below I c . In the second 
case the actual switching procedure is mainly noise in- 
duced. A suitable description of switching times in this 
regime can be obtained from the Fokker-Planck approach 
which was recently introduced by Apalkov and Visscher 
in the context of thermal fluctuations^^. Within this 
approach switching rates are specified by an Arrhenius 
like law with an effective temperature T e g. The latter 
differs from the real temperature T, as it is influenced 
by the damping and the spin torque. In the sequel we 
present a generalization of the method to nonequilibrium 
noise, and show that the spin shot noise alters the effec- 
tive temperature. 

Let us start our consideration with the Fokker-Planck 
equation as introduced by Browr*^. We denote the prob- 
ability density for the magnetization of a mono-domain 
particle by p(m,t). The corresponding Fokker-Planck 
equation can be written in the form of a continuity equa- 
tion 



dp(m,t) 

at 

with probability current 



-V-j(m,i) 



(32) 



Here rhdet denotes the deterministic part of the stochastic 
LLG (UJ and D is the random field correlator. We recall 
that the dynamics governed by (P) conserves the absolute 
value of m. As a consequence the movement of the tip of 
m is restricted to the surface of a sphere, which we will 
call the m-sphere. The gradient and the divergence in 
(1521 and (1331) are 2-dimensional objects, both living on 
the m-sphere. 

We now observe that in presence of anisotropy the 
phase space will be in general separated. The poten- 
tial landscape will exhibit different minima referring to 
stable and meta-stable states of the magnetization. Pre- 
cession of the magnetization takes place around one (or 
more) of these equilibrium positions. We refer to orbits of 
constant energy as Stoner-Wohlfarth (SW) orbits. Now, 
considering the dynamics of the magnetization vector one 
can distinguish two different time scales. The time scale 
for the angular movement, on the one hand, is charac- 
terized by the precession frequency. On the other hand 
there is also a time scale for a possible change in energy. 
In the following we will require that the time scale for the 
change in energy is much longer than the time scale for 
constant energy precession. In other words: we assume 
that the magnetization vector stays rather long on a SW 
orbit before changing to higher/lower energies. In this 
low damping and small current limit we can introduce 
an energy-dependent probability density by identifying 
p'j(-E(m),£) = p(m,t), where the index i takes into ac- 
count that the energy dependence may be different in 
different regions of the m-sphere. The above mentioned 
time scale separation allows us to average out the move- 
ment along the SW orbit and to be concerned with only 
the long time dynamics. 

The idea of the FP approach is to translate equation 
(f3"2"j) into a corresponding equation for p'AE). For ther- 
mal noise this has been done in Ref[14|. We now give 
a generalization of the method to the angle-dependent 
spin shot noise of section [TTJ To this end we write the 
correlator (j2"6")l in the form 



D{9)=D th + D [l-Pcost 



(34) 



where D t h — 7 ^ B v is the thermal part and Dq[1 — 
Pcosff] the nonequilibrium part of the correlator M' 

We abbreviated the angle-independent part of the spin 
shot noise by D - We also used 



P = 



Gp — Gap 
Gp + Gap 



(35) 



In general we can write the Fokker-Planck equation for 
the distribution p' i (E(m),t) in the forir*^ 



1 P l {E)dp' i {E,t) 



M s p 



Of 



= ~jf(E,t), (36) 



j(m, t) = p(m, t)m det (m) - DVp(m, t) 



(33) where Pi{E) is the period of the orbit with energy E. j\ 



is the probability current in energy. It is given by 
jf(E,t) = fl}(m,t) x dm] ■ m = - 1 ap' l (E,t)lf(E) 



■7Vi(£,i)m p -ir - 



dp'(E) 
dE 



s^thi-e^i 



M„D, h I, 



(37) 



The constant J is defined in such a way that Jm p = 
7/(MV)I s if m p is a unit vector in direction of I s . 
Furthermore we have introduced the following integrals 
along the SW orbit 



I 



D th 



If - P <b cos BHeffdm 



lf(E) = ^H eff dm 
lf(E,t) = l dmx m 



(38) 
(39) 

(40) 



A steady state solution of the FP equation is obtained 
by setting jf = 0. From (|3"Tf we get the following differ- 
ential equation for the probability density p\: 



7 



D th M, 



■Xiffli-a + THWJIs-VttiE) 



dE 

( 41 ) 
where the right hand side serves as a definition of an 

inverse effective temperature f3[{E). From (|4Tj) one can 
see that, depending on the sign of the spin current, the 
spin torque may either enhance or diminish the damp- 
ing, leading to a lower or higher effective temperature, 
respectively. In (|4*Tj) we have defined 



Vi(E) = 



m„ 



If(E) 



lf(E) 



and 



\{E) 






(42) 



(43) 



rji can be viewed of as the ratio of the work of the Slon- 
czewski torque to that of the damping^. The quantity A 
gives the renormalization of the effective temperature as 
compared to the pure thermal case. We can write for A 



HE) = ^ 

cff 



(44) 



where T c g is the effective temperature when only equilib- 
rium noise is present, and T c ' ff the effective temperature 
when both, equilibrium and nonequilibrium noise are in- 
cluded. It should be observed from (|4"Tj) that the effective 
temperature is in general energy dependent. The corre- 
sponding probability distribution will thus, in general, 
differ from the Boltzmann distribution. However, when 
we turn off the nonequilibrium, X(E) = 1 and J = 0. 
In this case the solution of (|4"Tj) is exactly a Boltzmann 
distribution. 
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FIG. 2: (Color Online) A = T e ff/T' e }f as a function of c in the 
case m p j\|, e z for eV — k B T (red), eV — 5fcsT (magenta), 
eV = Wk B T (blue), eV = 20k B T (green). 



In the remainder of this section we evaluate A for an 
exemplary system with easy axis and easy plane anisot- 
ropy. The easy axis is chosen to be the z axis and the 
easy plane is the y-z plane. The magnetization direc- 
tion of the fixed layer, m p , is taken to be anti-parallel 
to the z axis. Let us use the following convention for 
the spherical coordinates: m x — cos$, m y — sin $ sin y), 
m z = sin d cos (p. The SW condition defines the orbits of 
constant energy. For our system it reads 



E(M) _ 



\H K M s (me z f + \M 2 s (me x f 



(45) 



We abbreviate n\ = HkMs (characterizing the strength 
of easy- axis anisotropy), n^ — Mj (characterizing the 
strength of easy- plane anisotropy) and d = &-, being 
the ratio of easy-axis to easy-plane anisotropy, so that 
(taking the magnetic constant fio = 1) we can obtain 
from equation ((45]) the dimensionless energy c 



_ 2E 



= -dm. 



(46) 



This relation defines the 'potential landscape' of our sys- 
tem. We can distinguish three regions: Two potential 
wells, one around ip = (well 1) and one around ip = it 
(well 2), and a third region (region 3) with energies above 
the saddle point energy, separating the two wells. Switch- 
ing takes place if the magnetization vector changes from 
some orbit in the one well to an orbit in the other well. 
Equation (J4"6"l) defines the orbits of integration for the 
evaluation of (|43|. Let us concentrate on orbits lying in 
the potential well around ip = with energies c < 0^ 
In addition we assume a strong easy plane anisotropy, 
allowing to consider small deviations of i? around ^ . 

We fix the Gilbert damping to a = 0.01, the ratio of 
anisotropies to d = 0.028, the polarization to P = 0.81, 
and M s V/7 = 10H. These values define the ratio Do/ Dth 
as a function of eV/ksT. The results for A = T e ff/T^., 
are plotted in Fig. [2] As can be seen from the plot taking 
into account the nonequilibrium noise results in a renor- 
malization of the effective temperature. This renormal- 
ization is proportional to the applied voltage V and can 
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FIG. 3: (Color Online) A' as a function of c in the case 
m p U e z for P = 1 (red), P = 0.9 (green), P = 0.7 (blue), 
P = 0.5 (magenta), d = 0.028. 




be very strong for sufficiently large values of V. The de- 
viation from the purely thermal case (A = 1) approaches 
15% for eV = 5fcsT and is thus experimentally not negli- 
gible! For eV = 10&bT the deviation is even in the order 
of 25% and grows further with the voltage. The variation 
of A with energy is on the other hand very weak. This 
indicates that the influence of the angle dependence is 
rather small or in other words: The angle-dependence of 
the correlator does not lead to a significant variation of 
T e ff with precession orbit. 

Let us continue our discussion of the renormalized ef- 
fective temperature by considering the limit where the 
equilibrium part of the correlator is much smaller than 
its nonequilibrium part and thus may be neglected. In 
this case we define the following quantity of interest 



\'(E) = ^-X(E) . 



(47) 



One should note the difference between A and A'. From 
(|4"4"|) we see that A is the ratio of the effective tem- 
peratures T e ff and T' e tt for systems without and with 
nonequilibrium noise respectively. On the other, from 
the definition (j47|) it is clear that A' is a measure for 
the influence of the angle dependence of the correlator. 
The stronger A' deviates from A' = 1 the stronger is the 
influence of the angle-dependence. 

In Fig. [3] we plot A' for our model system for d = 0.028 
and different values of P. As one can see from Fig. [3] the 
largest deviation from X'(E) = 1 (corresponding to the 
strongest influence of the angle-dependence) is present 
at the minimum of the well (c = — d = —0.028). The 
smallest deviation from X'(E) — 1 is observed for orbits 
which lie near the separatrice. The overall change of 
\'{E) for P = 1 is of the order of 10%. 

These results provide a good insight into the influence 
of the angle- dependence. As A' ~ l/T'. », a small value of 
A' indicates a "hot" spot whereas large values of A' corre- 
spond to " cold" spots on the m-sphere. For the particular 
system under consideration, c/. Eq. (|45|) , the equilibrium 
position of the magnetization is roughly along the z axis. 



FIG. 4: (Color Online) Hot spots (red) and cold spots (blue) 
on the M-half-sphere in case of m p t4- e z ■ The noise intensity 
is highest at the bottom of the well. 



SW orbits of precession are symmetric with respect to 
this axis. At the bottom of the well 9 — tt and the spin 
shot noise has its maximal value. We hence expect a hot 
spot at the minimum of the well. With increasing en- 
ergy the orbits will become larger. The angle 9 will vary 
along these orbits. However as the orbit energy grows the 
trajectories increasingly go through regions of smaller 9, 
so that the average value of 9 will diminish with orbit 
energy. As a consequence the nonequilibrium noise will 
become smaller as well. Cold orbits should be hence those 
that are in the vicinity of the separatrice. This is exactly 
what can be read off from Fig. [3J Our findings are thus 
in agreement with the geometrical situation. Cold spots 
and hot spots on the m-sphere are shown in Fig. 01 



IV. SWITCHING TIME OF SPIN-TORQUE 
STRUCTURES 



The switching process can be analyzed by performing 
numerical simulations of the Langevin equations of mo- 
tion with the inclusion of temperature and shot noise 
via the random field term. In this section we present 
such simulations for Gilbert damping of a = 0.01, an 
anisotropy ratio of d = 0.028, and with a spin torque cur- 
rent characterized by J and polarized in the m p = — e z 
direction. 

Before going further though, it would be useful to con- 
sider how the system acts in the absence of the noise. In 
such a case the switching occurs when the energy cur- 
rent (1.34) is positive for all values of energy between 
the starting position (say positive z direction) and the 
saddle point. Since the probability function pf(E.t) is 
always positive it stand to reason that a switch will only 



happen if < -aIf{E) 
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FIG. 5: (Color Online) Plots A(E) = -^lf(E) + -^-m p ■ 

if 1 (E) as a function of energy c = 2E/K2 over a range of spin 
torque current. From bottom to top: Blue = 0.77 J c , Red = 
1.08 J c , Yellow = 1.55X, Green = 3.10J C 



quantity as a function of energy for various values of the 
spin-current J in Fig. O From this we also gain a useful 
reference value for the critical current current which is 



J, 



aI?(E 3ad ) 



,j ,. — 0.00645M S . The positive value signi- 

fi.es the tendency towards the switching. In the first ex- 
ample with J = 0.77 J c the noiseless system, being driven 
by the dissipation towards the stable position, does not 
switch. It is worth noticing that in the presence of the 
noise the switching nevertheless does occur, but it takes 
exponentially long time. In the three other examples 
J > J c and the magnetization current is always directed 
towards the saddle. Therefore even the noiseless system 
does switch and the noise serves to introduce an uncer- 
tainty in the switching time. 

Putting the thermal noise back into the system, we 
set the noise strength parameter to D t h = 0.0000l7M s . 
Simulations are then run by starting each particle at 9 = 
0, allowing it to come into thermal equilibrium with the 
system, turning the current on and calculating how long 
it takes for it to go past the saddle point into the second 
well. This is done for many particles for a given current 
value and over several different current values. A typical 
trajectory of the system is represented by the graph 9 as a 
function of time in Fig. [6]for J = 1.08 J c . It may be seen 
that it takes many revolutions before the system finally 
switches to the basin of attraction of the true stationary 
points at t sw 700. 

Let us estimate the contribution of the nonequilibrium 
noise to the switching process under realistic experimen- 
tal conditions. From Eq. (|26|) we obtain the relationship 
between equilibrium noise and non-equilibrium noise 
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FIG. 6: A typical realization of 
units of (-yMs)^ 1 ) for J = 1.08J C . 



as a function of time (in 
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FIG. 7: (Color Online) (Left) Shows the switching probability 
as a function of time (in units of {■yM a )~ 1 ) for various current 
values. (From left to right: Green = 3.10J C , Yellow = 1.55J C , 
Red = 1.08J C , Blue = 0.77J C ). (Right) Shows the average 
switching time (in units of (7M1,) ) as a function of 4-. 



where J c is the minimum spin current needed to cause a 
switch in the absense of noise, we obtain 



A, 



D 



CC| 



h-y J c 
2k B a T 



(50) 



Using the material parameters in Ref. Q, M s — 1440 
emu, nanopillar volume = 1.97 x 10~ 17 cm 3 , and switch- 
ing current I c « 10 9 A/ cm 2 , the shot noise, Do, at the 
critical current is equivalent to the thermal noise, D t h, 
of temperature T c ps 15 K. Even though the theoreti- 
cal switching current and the experimentally determined 
switching current of Ref. [6| differ by more than an order 
of magnitude (such a discrepancy is mentioned in many 
experimental studies, including Ref. [6j) the theoretical 
approach can be used to gain insight into how the shot 
noise at the switching current scales with the parameters 
of a nanopillar. 

If we set Dq = D t h again and include our calculation 
of the switching current we find T c oc J c /{aP) where 
J c is the critical spin-torque current needed to cause a 

(48) spin-flipping event (in the absence of noise) and P is the 
degree of polarization of the current. J c is determined 
by the relative values of Hk and M s and has the units 
of magnetization. For the simple case with only easy- 
axis anisotropy (E — —^fi^H^Mg cos 2 9), J c — aH^, and 

(49) therefore T c ex -#. Therefore for materials and/or geom- 



etry with a bigger anisotropy field the shot noise may be 
a dominant source of noise at temperatures well above 
15K. Moreover, in the simplest case where we only have 
uniaxial anisotropy and an applied external held along 
the easy axis it can be shown that 



J c = a{H k + H c 



(51) 



This means the relative strength of the non-equilibrium 
noise scales with the potential well depth of our system. 
By applying an external held we can increase the maxi- 
mum allowed current before a switch takes place and thus 
increase the importance of non-equilibrium noise. 

Since the initial condition is taken out of a stationary 
distribution (without the spin current) and subsequent 
evolution is subject to the Langevin noise, the time of 
the switching is a random quantity. The percentage of 
trial systems that have switched as a function of time is 
shown in the left panel of Fig. [7] for four different values 
of the spin-current. The time derivatives of these graphs 
provide probability distribution functions of the switch- 
ing time. One may then evaluate the first moment of 
these distributions which gives the mean switching time 
for a given value of the spin current. The right panel of 
Fig. [7] shows such a mean switching time as a function of 
J/J c . One may notice that for J < J c the switching time 
grows exponentially, while for J > J c the switching time 
becomes relatively short (although it is still substantially 
longer than the inverse precession frequency). 



V. CONCLUSION 

We would like to conclude with the following remarks. 
The study of noise in dynamical magnetic systems is a 
broad and fascinating held. In particular, in view of po- 
tential applications of magnetic nanodevices, both equi- 
librium and nonequilibrium noise may play an impor- 
tant role. For example, the stability of magnetic stor- 
age devices is strongly influenced by thermal fluctua- 
tions. The functionality of new generation technologies 
(such as the magnetic random access memory (MRAM) 33 
with STT writing or the racetrack memory 3 ^) is largely 
based on the spin torque phenomenon. The latter is a 
nonequilibrium effect and thus besides the temperature 
also nonequilibrium sources of noise may play an impor- 
tant role. 



A way to introduce fluctuations into magnetization dy- 
namics is to add a random component to the effective 
field or to the current in the phenomenological LLG equa- 
tion. The noise is then defined by the value of its correla- 
tor (and its higher order cumulants). The determination 
of the noise correlator is of great importance as it defines 
the noise properties. In addition it may give insight into 
the physical context. 

A powerful and very flexible tool is the Langevin ap- 
proach based on the Keldysh path integral formalism. 
Starting from a microscopic model one derives the equa- 
tions of motion for the magnetic system. Fluctuations 
naturally arise as a generic feature of the Keldysh ap- 
proach. We have demonstrated the applicability of this 
method to magnetic systems on the example of spin shot 
noise in magnetic tunnel junctions. The spin shot noise 
correlator arose naturally, as a consequence of the sequen- 
tial tunneling approximation, in second order in spin flip 
processes. 

The Keldysh formalism is however not restricted to the 
system described above. In particular it may be used in 
the context of nonuniform magnetic textures (as domain 
walls for instance). Promising advances in this direction 
have already been reported 3 ^ and demonstrate the ver- 
satility of the method as well as inspire us with curiosity 
about future developments. 

Finally, to investigate the influence of the spin shot 
noise on spin torque switching rates we have generalized 
the Fokker-Planck approach of Ref 14| . We have shown 
that the nonequilibrium noise manifests itself in a renor- 
malized effective temperature. In particular at low tem- 
peratures we could observe a significant variation of the 
noise with orbit energy, reflecting "cold" and "hot" tra- 
jectories of the magnetization vector with respect to the 
noise intensity. 
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